extern doublereal dzsum1_(integer *n, doublecomplex *cx, integer *incx);
